Generating Schrodinger-cat states in momentum and internal-state space from 
Bose-Einstein condensates with repulsive interactions 



CO 

o 
o 

CN 

» 

o 
O 

CO 
CN 



43 

o 



CO 



I 

C 

O 

o 



> 

CO 

in 
o 

CO 

o 

c3 



o 
o 



X 
S3 



J. Higbie and D. M. Stamper-Kurn 

Department of Physics, University of California, Berkeley CA 94720 
(Dated: February 2, 2008) 

Resonant Raman coupling between internal levels induced by continuous illumination of 
non-collinear laser beams can create double-well momentum-space potentials for multi-level 
"periodically-dressed" atoms. We develop an approximate many-body formalism for a weakly inter- 
acting, trapped periodically-dressed Bose gas which illustrates how a tunable exchange interaction 
yields correlated many-body ground states. In contrast to the case of a position-space double well, 
the ground state of stable periodically-dressed Bose gases with repulsive interactions tends toward 
a Schrodinger cat state in the regime where interactions dominate the momentum-space tunnelling 
induced by the external trapping potential. The dependence of the momentum-space tunnelling 
and exchange interaction on experimental parameters is derived. We discuss how real-time control 
of experimental parameters can be used to create Schrodinger cat states either between momentum 
or internal states, and how these states could be dynamically controlled towards highly sensitive 
interferometry and frequency metrology. 

PACS numbers: 03.75.Gg,05.30. Jp,52.38.Bv 



Following our mastery over the internal and external 
states of individual atoms, the scientific frontier advances 
to the full control over the quantum states of many-body 
systems. Entangled states, in which the constituents of 
a many-body system display non-classical correlations, 
play a key role in our developing understanding of quan- 
tum information and decoherence, and may find practical 
use in quantum communicationlll, quantum computing, 
and high-precision metrology S S Various tech- 
niques are being developed that generate entanglement 
deterministically in, for example, spin-squeezed atomic 
ensembles 0, trapped ions 0, the electromagnetic field 
@, and superconducting circuits [^.Hof. 

Ultracold neutral atoms offer a promising route to gen- 
erating highly entangled many body states. Schemes 
for generating such entanglement rely on interatomic in- 
teractions which provide "non-linear" elements as seen 
from the viewpoint of single-particle dynamics. Such 
non-linear terms are provided naturally through binary 
collisions between ground-state atoms, as utilized con- 
trollably by Mandel et al. [Tlj, and can be further en- 
hanced through molecular resonances or by the assistance 
of near- resonant photons |12l 1 1 3| . 

A paradigmatic system in which interatomic interac- 
tions induce entanglement is a collection of ultracold in- 
teracting bosons in a parity-symmetric double-well po- 
tential El H E E3 • The lowest -energy single-particle 
states are the parity-even ground state \ipo) and parity- 
odd first excited state \tpi), separated in energy by the 
tunnel splitting J. If tunnelling and interactions are 
weak with respect to the energy spacings to other single- 
particle states, the many-body system may be described 
by the two-mode Hamiltonian 

~ f4cL + c[c B ) - UN R N L (1) 
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whence terms dependent only on the total number of 
atoms in the system are omitted. Here cr and cr, denote 



annihilation operators for particles in the right- = 
[|Vo> + IV>i}]/\/2) or left- (\L) = [|Vo> - \*l>i)]/V5) well 
states, respectively. The number operators Nr = CrCr 
and Nl = c^Cl count particles in the right or left wells, 
respectively, of the double-well system. The parameter 
U gives the energy due to short-range interactions of a 
pair of atoms located in the same well. 

This simple Hamiltonian leads to highly-correlated 
many-body ground states through the interplay of tun- 
nelling and interaction. One finds three limiting behav- 
iors. If the tunneling rate dominates, the many-body 
ground-state of N bosons is driven to the factorized, 
uncorrelated state in which all N atoms are identically 
in the single-particle ground state \ipo)- In the limit 
N\U\/J 3> 1, the interaction energy dominates. For 
repulsive interactions (U > 0), a many-body state di- 
vides itself evenly between the two potential minima, 
generating the state |^) cx {$ R ) N / 2 ($ L ) N / 2 \Q) , where |0) 
is the vacuum state. Experimental evidence for such 
a "number-squeezed" state has been obtained by Orzel 
et al. in a many- well potential |18| . Similar physics is 
responsible for the observed Mott-insulator phase in a 
three-dimensional periodic potential |19| . 

For attractive interactions (U < 0), the ground-state is 
quite spectacular: a Schrodinger-cat state \$?) cx [(c R ) N + 
(4)^10) formed as a superposition of states in which all 
atoms are found in one of the two wells. The genera- 
tion of Schrodinger cat states is a tantalizing goal; at 
present, the largest atomic Schrodinger-cat states con- 
tain just four atoms [?J. However, accomplishing this 
goal by imposing a double-well potential on a collection 
of bosons with attractive interactions is a daunting task: 
such Bose-Einstein condensates are unstable to collapse 
20] , implying that the number of bosons placed into the 
superposition state will be small; moreover, exceptional 
spatial control over trapping potentials and over the in- 
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teraction strength is required. 

The above discussion pertains to atoms for which the 
position-space potential is a double well, while the kinetic 
energy term p 2 /2m may be regarded as a momentum- 
space harmonic "potential. " In this work, we point out 
that an analogous many-body system can be crafted, in 
which the roles of position and momentum are inter- 
changed. As discussed in previous work 0, 152, l23| . 
the dispersion relation of multi-level atoms placed in a 
spatially-periodic coupling between internal states can 
take the form of a momentum- space double-well "po- 
tential. " Addition of an harmonic position-space trap- 
ping potential produces a double-well system dual to 
that discussed above. Assessing the roles of interparti- 
cle interactions in this situation, we find that the many- 
body ground states for bosonic atoms in such a system 
are highly- entangled in both momentum- and internal- 
state space. In opposition to the situation of bosons in 
position-space double wells, maximally-entangled states 
may be generated in the case of repulsive interactions, al- 
lowing the creation of such states starting with large, sta- 
ble Bose-Einstein condensates while obviating the need 
for exacting control over potentials with extremely small 
spatial dimensions. 

The roles of repulsive or attractive interactions in 
determining the behaviour of interacting bosons in 
momentum-space double-well potentials are reversed 
from those in the case of position-space potentials. As 
discussed below, this reversal is due to an exchange term 
which arises in the evaluation of the interaction energy 
of a Bose gas that occupies several distinct momentum 
states. For repulsive interactions, for example, this ex- 
change term disfavors a macroscopic occupation of more 
than one single-particle wavefunction (the so-called frac- 
tionation of a condensate as discussed by Leggett 
and others), thus favoring a Schrodinger-cat superposi- 
tion. We find that the strength of the exchange term 
can be dynamically varied by varying parameters of the 
periodic coupling field which generates the double-well 
potential. 

Following a derivation of the momentum-space double- 
well potential (SecQJ, we develop a two- mode approxi- 
mate treatment of this system which provides expressions 
for the interaction and tunnelling energies and thereby 
clarifies their dependence on experimental parameters 
(Sec. [njl. Similar work by Arecchi and Montina [2^ 
treats this system using the Gross-Pitaevskii equation as 
the starting point, a numerical scheme which identifies 
the onset of correlated many-body ground states. Their 
work appears to reproduce our analytic approach in the 
limit of weak interactions, while the validity of their ap- 
proach (or, indeed, any extant treatment) for stronger 
interactions is difficult to establish. Furthermore, in Sec- 
tion lllll we show that maximal entanglement can be gen- 
erated purely between internal states or between momen- 
tum states, offering a route to Heisenberg-limited atomic 
clocks 3] or atomic interferometry 0, respectively. 



I. ORIGIN OF THE MOMENTUM-SPACE 
POTENTIAL 

We consider bosonic atoms of mass m with two inter- 
nal states, | A) and \B), at energies Hu>a and hu>B, respec- 
tively. These atoms are exposed to laser fields of frequen- 
cies ui\ and u>2 = — 0J2) an d wavevectors ki and k2 
which induce Raman transitions between the two inter- 
nal states (see Fig. QJi,b). Thus, a Raman process from 
state I A) to \B) imparts a momentum of Kk = h(k.i — k.2) 
and kinetic energy of h5 = fou> — h(ujB — ^a)- The con- 
tinuous Raman coupling can be regarded as a spatially- 
periodic coupling field between the internal states of the 
atoms, represented by an off-diagonal potential Vr — 
Ml/2 (e-^ k ' r - w *)|A)(S| +e i( - k - r - u) V\B){A\). Here fl is 
the two-photon Rabi frequency, taken to be real, which is 
determined by dipole matrix elements, the detuning from 
intermediate resonances, and by the laser intensities. 

It is convenient to analyze this constantly-driven sys- 
tem in terms of "periodically-dressed" states, which are 
coherent superpositions of both internal and external 
(momentum) states. As developed in Refs. 0,H2, this 
treatment reveals marked anisoptropy and tunability in 
the superfluid characteristics of a Bose-Einstein conden- 
sate formed of such a dressed-state gas. Given a two- 
component spinor wavefunction ^(r) = (^(r), ^s(r)) 
to describe single-atom states in the {\A), \B}} basis, we 
transform to a frame which is co-rotating and co-moving 
with the driving laser fields, yielding a spinor wavefunc- 
tion ^(r) with components iPa,b(*) — exp(±i(k ■ r — 
Lot) /2)V ) A,s(r). This transformation yields a Schrodingcr 
equation with two deBroglie wave solutions at each 
momentum ftq, which we call the the periodically- 
dressed states. The two-branch dispersion relation of the 
periodically-dressed atoms has the form (using 7r for plus 
and mu for minus) 

fiu* tlt (q) = q 2 + \ ± \ yWk-^ + fF (2) 

where energies and wavevectors are scaled by the Raman 
recoil energy Ej~ = h 2 k 2 /2m and wavevector k, respec- 
tively. 

Under the conditions of exact Raman resonance (8 = 
0) and sufficiently small Rabi frequency f2, a degeneracy 
of momentum ground states occurs as the lower disper- 
sion relation takes the form of a double-well potential, 
with minima at ±Q/2 (Q ~ k) (Fig.[TJ;). This potential 
can be quickly tuned by modifying the laser based pa- 
rameters. The detuning from Raman resonance 8 breaks 
the ground-state degeneracy, favoring the right- or left- 
well states. The spacing between the potential wells is 
controlled by the Raman momentum transfer Hk, which 
can be varied by reorienting the laser beams. Finally, the 
Rabi frequency £1 changes both the height of the barrier 
between the wells, as well as the internal-state character 
of states on either side of the well [z| • 

In the case of a position-space double well, the kinetic 
energy of the atoms forbids a complete localization of 
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FIG. 1: Creating a double-well momentum-space potential 
with Raman excitation, (a) Laser beams of frequency u)i and 
u>2 may induce Raman transitions between internal states \A) 
and \B). 5 = {u>\ — uj?) — ojq is the detuning from the Raman 
resonance, (b) Such a Raman transition imparts a momentum 
transfer of ftk. = ft(ki — IC2), where ki and k2 are the wavevec- 
tors of the Raman coupling lasers, (c) Atoms exposed to con- 
tinuous Raman excitation can be described by a two-branch 
dispersion relation. Energies (scaled by Ek) for the lower 
{hujp) and upper (fia) T ) dressed states are shown for wavevec- 
tors (scaled by k) in the direction of the Raman momentum 
transfer. In the case of exact Raman resonance (S = 0) and 
small Raman Rabi frequency (shown for htt/Ek = 1/8), the 
lower dispersion relation takes the form of a parity-symmetry 
double-well potential. 



atoms in either of the wells, thus introducing tunnelling. 
Similarly, in our case of the momentum-space double- 
well, a tunnelling between well-defined momentum states 
can be induced by adding a spatially-dependent term 
to the Hamiltonian. We consider adding an internal- 
state independent position-space trapping potential of 
the form V(r) — mw ( 2 r 2 /2 |23 |. The Hamiltonian for this 
system may be written in the basis of the periodically- 
dressed states introduced earlier, using the expansion 



V>(r, t) 



rf 3 q 

(2tt) 3 / 2 



R(q) 



7r(q) 

M(q) 



(3) 



where R(q) = e lcr v e (i)/ 2 where a y is a Pauli matrix and 
the mixing angle 0(q) is denned by the relation cot 0(q) = 
(5 — fiq- k/m)/fi. The wavefunction 0(q) = (7r(q), //(q)) 
in the space of periodically-dressed momentum eigen- 
states obeys a Schrodinger equation with Hamiltonian 



W+(q) 
w_(q) 



1 



^Rt(q)V^R(q) (4) 



Hpd = ft 

The position-space trapping potential is seen in mo- 



mentum space as a kinetic-energy-like term (involving 
Vg), which accounts also for the variations with q of 
the periodically-dressed eigenstate basis. The relevance 
of this kinetic-energy-like term is measured by the di- 
mensionless effective mass parameter M, which is re- 
lated to the Lamb-Dicke parameter 77 = y/hk 2 /2rnu>t 
as M = 2rj 2 . For weak spatial confinement (M ^> 1), 
the low-energy single-particle states are restricted pri- 
marily to the lower periodically-dressed states, with the 
two lowest states split by a small tunnelling energy J. 
For strong confinement (M <C 1), the single-particle 
states are admixtures of upper- and lower-dispersion- 
rclation periodically-dressed states, and a simple double- 
well treatment is no longer adequate. 

To introduce interatomic interactions, we make use of 
the field operator 0(q) = (n(q), /t(q)) the components 
of which annihilate particles in the lower (/1) or upper 
(tt) periodically-dressed states. These operators obey 
Bose commutation relations. We consider low-energy, 
binary, elastic collisions with a state-independent scat- 
tering length a. The interaction Hamiltonian then takes 
the form H mt = f (/ d 3 qn(q)n(—q}^j , neglecting terms 
dependent only on the total number of atoms N, where 
n(q) is the spatial Fourier transform of the density op- 
erator, and g — 8irka is the properly scaled interaction 
parameter. The density operator is given as 0, 12^ 



n(q) 



d 3 K 

(2^ 



3)Rt( K + |) R(K _|^( K _3) 

(5) 



II. THE TWO MODE APPROXIMATION 

To simplify our treatment, let us consider only the sit- 
uation in which the two lowest energy eigenvalues (in 
the absence of interactions) are well separated from the 
remaining energies, and thus a two-mode description of 
the many-body system is adequate. The two lowest- 
energy single-particle states have wavefunctions </>o(q) 
and 0i (q), mode operators Cq and &x, and energy split- 
ting J. We define the right- and left-well states as 

0R,i(q) = (<Mq) ± <?i(q)) A/2, and the mode opera- 
tors as cr^l = (co ± di)/y/2. Under this approximation, 
we can evaluate the density operator n(q) by expanding 
the field operators fi(q) and 7r(q) in the basis of energy 
eigenstates, and then truncating the expansion after the 
first two states. We then express the density operator as 
«(q) - Ei,,--A4;(q)cjcj- with 



A^(q) l 



d 3 K 

(2tt) 3 / 2 



$( K +|)R(k+|)R( K -|)^( K +|) 

(6) 

with indices i,j € {R, L} denoting cither the right or 
left states. Accounting for prop erties of the right- and 
left-well states under parity 26], we obtain 

n(q) ~ Af RR ((i)N + Ar RL (<l)c f R c L + A/ki(-q)4ci? (7) 
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where N is the total number of atoms in the system. 
Dropping terms that depend only on N, one thus finds 



H 



hlt ~ -UN R N L + ^(c R c L +c\c R 



' R c L c L + & l & l crCr 



+j 2 (44 

where the energies U, J\ and Ji are given as 



U 



RL 



Jl = 2Ng J d 3 qM RR (q)Af RL (q) 
J2 = § / d 3 qA^ L (q)A^ L (-q) 



(8) 

(9) 
(10) 

(11) 



Before delving further into the implications of the in- 
teraction Hamiltonian given above, let us consider the 
weak confinement (large M) limit at which the right- 
and left-well states contain no population in the upper 
(71") dressed states, are Gaussian functions well-localized 
(in q space) at the right and left potential minima at 
±Q/2, with rms widths a q ~ (2M)~ X / 2 , and have av- 
erage density per particle (n) = a^/ir 3 / 2 . In calculating 
the integrals of Eq. we may thus use local values of 
the rotation matrices R(q) at ±Q/2, as appropriate. We 
then find 



U 
J-2 



-g{n) sin 



2Ng(n)sm6 (^j x e ^ 
|(n> sin 2 0(^1 < 



(12) 

(13) 
(14) 



These approximations can be further simplified by setting 
Q = k, and thus sin0(Q/2)) ~ ^^/(l + fl 2 ). 

Let us presume that the interaction-mediated tun- 
nelling terms J% and J2 are negligibly small. We thus 
recover the many-body Hamiltonian of Eq. ^ which had 
applied to the position- space double well. However, exam- 
ining the interaction energy U, we now find that the roles 
of attractive and repulsive interactions are interchanged 
between the position-space and momentum-space double 
well treatments: for the case of the momentum-space 
double-well potential, repulsive interactions (g > 0) 
yield U < and thus are the source of Schrodinger- 
cat ground states, while attractive interactions (g < 0) 
yield U > and thus are the source of number-squeezed 
ground states. This opens the door to the production 
of Schrodinger-cat states for large, robust Bose-Einstein 
condensates with repulsive interactions through the use 
of binary collisions as a non-linear coupling. 

This reversal of roles for repulsive and attractive in- 
teractions in momentum space can be understood in the 
context of a single-component Bose-Einstein condensate 
in a harmonic potential. Repulsive interactions lead the 



ground-state condensate wavefunction to become larger 
spatially, thus causing the momentum-space wavefunc- 
tion to contract - appearing as an attractive interaction 
in momentum-space. Equivalently, one may note that 
in a Hartree approximation to the many-body state of a 
Bose-Einstein condensate, the interaction energy is eval- 
uated as an interaction energy density proportional to 
the square of the density. The density of a Bose gas 
occupying two distinct momentum states would be spa- 
tially modulated by the interference between the momen- 
tum states, and, therefore, its interaction energy oc gn 2 
would be increased (decreased) for the case of repulsive 
(attractive) interactions. Thus we find repulsive interac- 
tions leading to Schrodinger-cat states, a superposition 
of states in which only one single-particle state is macro- 
scopically occupied, while attractive interactions would 
lead to a "fractionation" of the condensate between dis- 
tinct momentum states. 

In the system we have described here, the interaction 
energy U is a tunable exchange term which, in the case 
of repulsive interactions, suppresses the superposition of 
atoms in two distinct momentum states and thereby fa- 
vors the maximally-entangled state. The exchange term 
in the interaction Hamiltonian arises from the presence 
of atoms in identical internal states, but different mo- 
mentum states. By adjusting the strength of the Rabi 
coupling f2, the admixture of internal states \A) and \B) 
in the right and left well states is varied. For small SI, 
atoms in the right- well state are almost purely in the \A) 
state, and atoms in the left-well state are almost purely 
in the \B) state. Thus, the strength of the exchange 
term is suppressed (sin 2 #(Q/2) <C 1), and the product 
state I*) cx (c R + cl) n \0) remains energetically favored. 
For larger fi, the right and left-well states both contain 
admixtures of the two internal states; e.g. the right-well 
state contains atoms in state \A) nearly at rest, while the 
left-well state contains atoms in state \A) at momentum 
~ — ftk. Significant exchange terms now suppress the 
product state in favor of the maximally-entangled state. 

The dependence of the many-body ground state on the 
ratio U/J and the atom number N has been worked 
out by various authors E E3 for the position- 
space double-well potential. As we have obtained a 
Hamiltonian identical in form to the position-space treat- 
ment, these predictions would apply directly to our 
scheme. Significant deviations from the factorized, non- 
interacting atom solution begin at \U/J\ > 1/N; for 
repulsive (U > 0) or attractive (U < 0) interactions 
this would shift population either away from or to- 
wards the state \N L = N/2,N R = N/2), respectively. 
Perturbation analysis shows that the many-body state 
will be well approximated by the Schrodinger-cat state 
(|JV,0) + \0,N))/V2 when U/J > 1/VN. 

Now returning to Eq. we see that two extra terms 
appear in Ti.- m ^ which do not conserve the number of par- 
ticles in the right- and left- well states, representing a form 
of interaction- mediated tunnelling. The first, involving 
jW, modifies the tunnelling energy between the right 
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and left wells. The effect of this term is to reduce the tun- 
nelling rate for repulsive interactions, and increase it for 
attractive interactions. This can be understood by noting 
that repulsive interactions tend to raise the energy of the 
single-particle ground state more than the first excited 
state since the ground state has a higher density. Our 
previous analysis of the two-mode approximation accom- 
modates this term by a redefinition of J as J — > J — 
leading to no further complications. The second tun- 
nelling term, involving J^ 2 \ describes collisions which 
can redistribute two atoms from the right well to the left 
well, and vice versa. While this term modifies the con- 
clusions of our simple treatment, we have seen that the 
magnitudes of both JW and are exponentially sup- 
pressed with respect to U, and thus the dominant role of 
interactions is to create the aforementioned many-body 
ground states. These terms will, however, play an impor- 
tant role in determining the strength of interactions at 
which many-body correlations begin to become evident; 
the scaling of J (independent of g), and U with 

the control parameters (f2, g, M, etc.) are different and 
thereby provide more tunability to the system. 

We performed numerical calculations of the lowest- 
energy single-particle states to verify the simple scal- 
ing behaviour of various terms in the Hamiltonian de- 
scribed above. The single-particle Hamiltonian of Eq. 
0] is separable in the three Cartesian coordinates de- 
fined so that one axis (say z) lies along the direction 
of the Raman momentum transfer k. We then have 
the lowest energy eigenstates as </>o,i(q) = 4>o,i(Qz) x 

exp(— (q% + qy)/4(Tq) /^J2na 2 : a product of harmonic- 
oscillator ground states in the x and y directions and nor- 
malized solutions 4>o,i(q z ) to the one-dimensional, two- 
component Schrodinger equation derived from Eq. 

These one-dimensional eigenstates, calculated using a 
restricted basis set of Fourier components over the do- 
main — 3 < q < 3, are shown in Figure [21 for two differ- 
ent values of the mass parameter M, and for the con- 
dition of exact Raman resonance 3 = and Rabi fre- 
quency Ml/Ek = 1/8. For weak spatial confinement 
(large M = 50), the lowest two energy eigenstates are 
indeed nearly entirely composed of lower dressed states 
(ir(q z ) ~ 0) and are well-approximated by the sum or 
difference of Gaussian right- and left-well wavefunctions 
centered at the potential minima near q z /k ~ ±1/2. For 
stronger spatial confinement (smaller M = 8), the en- 
hanced "momentum-space tunnelling" causes the wave- 
functions to become less confined in the potential min- 
ima. Concomitantly, a significant population appears in 
the upper dressed states. 

Parameters which enter into the two-mode Hamlito- 
nian are derived from the numerically-obtained wave- 
functions, and are shown in Figure One finds the scal- 
ing J ~ exp(— M/4) as one expects [27|]. The numerical 
results confirm the scaling behaviour for found in 
the weak confinement (large M) limit fEas. 1121— 114|) . but 
differ slightly from the scaling predicted for since 
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FIG. 2: The two lowest-energy eigenfunctions (solid line for 
the ground state, dashed for the excited state) which define 
the right- and left-well modes for the two-mode approxima- 
tion. Component wavefunctions in the lower (fi(q z )) and up- 
per (ir{qz), shown in insets) periodically-dressed states are 
shown for M — 8 and M — 50. As the spatial confinement 
is weakened (larger M), the wavefunctions become further lo- 
calized in the minima of the double-well potential, and the 
population in the upper dressed state diminishes. Note the 
different parity of the fi and it components, as discussed in 
the text. Here hQ/Ek = 1/8, wavevectors are scaled by k, 
and the one-dimensional wavefunctions 4>o,i(q z ) are shown. 



the assumption that the right-well state remains Gaus- 
sian in the vicinity of the left-well potential minimum is 
incorrect. Nevertheless, is strongly suppressed for 
large M, approximately as exp(— M/3). Also shown is 
the energy splitting between the second and third lowest 
eigenenergies; one sees that for weak spatial confinement 
the lowest two energy eigenstates become well separated 
from the remaining eigenspectrum, establishing the va- 
lidity of the two-mode approximation. 

Finally, we stress that the two-mode treatment pre- 
sented here becomes invalid in the Thomas-Fermi regime, 
where the strength of atomic interactions dominates the 
zero-point energy in the confining potential. In this situa- 
tion, additional single-particle states must be considered, 
resulting in a complicated, self-consistent definition of the 
right- and left-well states which depends on the number 
of atoms in these states. How to properly treat such a 
situation remains an open question, and therefore an im- 
portant subject for experimental investigation. In related 
work, Montina and Arecchi [2j| obtain right- and left- well 
states through use of the Gross-Pitaevskii equation, but 
the validity of their treatment is questionable for systems 
which are not completely interaction-dominated. 
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FIG. 3: Numerical calculations of tunnelling and interaction 
strengths, a) The tunnel splitting J (open circles, solid line) 
becomes much smaller than the spacing between the second 
and third excited states (open triangles, solid line) for mod- 
erate values of M, establishing the validity of the two-mode 
approximation for weak confinement, b) Tunnelling energies 
J (open circles, solid line), j' 1 ' (open squares, dashed line) 
and J^ 2 ' (filled circles, dashed line) are suppressed for weaker 
confinement (exponentially with large M), while the strength 
of the momentum exchange energy U (X's, solid line) remains 
large. This provides a route to creating correlated many- 
body states adiabatically from uncorrelated states by gradu- 
ally weakening the spatial confinement. The energies in (a) 
are scaled by Ek- Dimensionless interaction parameters are 
plotted as -U/g(n), J (1) /2Ng(n), J (2) /(g{n)/2). A Rabi fre- 
quency liQ/Ek = 1/8 is chosen. 



III. APPLICATION TO HEISENBERG-LIMITED 
MEASUREMENT 

Several theoretical works have pointed out the poten- 
tial for correlated many-body states for improving the 
precision of atomic clocks Q and interferometers [faH^], 
in which phase shifts are measured at the Heisenberg 
limit A(f> oc 1/N, rather than the standard quantum limit 
Acj) oc 1/ VN, where N is the number of particles used in 
a single run of the experiment. In particular, Bollinger et 
al. described the use of the maximally-entangled state in 
such a measurement . In this section, we point out how 
dynamical control over experimental parameters can be 
used fruitfully to generate Schrodinger-cat states suited 
for the implementation of both Heisenberg-limited spec- 
troscopy and interferometry. 

In the system we have described, repulsive interatomic 
interactions lead to lowest energy many-body states 



which are the even- and odd-superposition Schrodinger 
cat states, \S±) = (|#r) ± |* L ))/\/2. The states 
and I^l) are distinguishable many-body states associ- 
ated with the right and left potential wells, respectively, 
which tend toward the limiting cases \*S?r) — ► |0, N) and 
I^l) — ► \N,0) for strong interactions. The energy sep- 
aration between these two states depends on the resid- 
ual overlap between \^r) and and the number of 
atoms in the system (derived, e.g., in Ref. ^H). In the 
weak-confinement limit M ^> 1, these right- and left-well 
states are comprised of the lower dressed states, which 
are superpositions in both internal and external degrees 
of freedom. Such a superposition appears to complicate 
the use of these Schrodinger-cat states for measurement 
applications. 

However, after these Schrodinger-cat states are formed, 
the dressed right- and left-well states can be converted 
to pure internal or external states adiabatically, i.e. by 
modification of experimental parameters on a time-scale 
which is fast with respect to tunnelling times so as to pre- 
serve entanglement, but slow with respect to timescales 
relevant to the Raman coupling (l/f2). If the Rabi 
frequency O is adiabatically lowered to zero (the Ra- 
man beams are slowly extinguished), atoms in the right- 
well state would be adiabatically converted to stationary 
trapped atoms in the \A) internal state, while atoms in 
the left well would be converted to stationary trapped 
atoms in the \B) state. Thus would be prepared a state 
ideal for Heisenberg-limited measurement of the internal 
state energy difference hujQ, potentially on a useful hyper- 
fine clock transition. Alternately, one could adiabatically 
ramp the detuning from the Raman resonance to large 
positive or negative values. For example, an adiabatic 
downward sweep of the frequency difference between the 
Raman laser beams (5 — > — oo) lowers the energy of the 
| A) internal state with respect to the \B) internal state in 
the rotating frame. This adiabatically converts the right - 
and left-well states to distinguishable momentum states 
with identical internal states (in this case the \A) state). 
One thereby prepares a state ideal for Heisenberg-limited 
atomic interferometry. 

Furthermore, the timed application of Raman coupling 
in this system can be used as a "magic beamsplitter" [2!| 
to prepare the Schrodinger-cat state dynamically. For 
simplicity let us make the identification \^r) = \0,N) 
and \^l) = \N, 0). Consider a system prepared in the 
\ty) = \N,0) state with no Raman coupling (f2 = 0) - this 
state corresponds to a zero-temperature Bose-Einstein 
condensate of N atoms in the \B) internal state. The 
Raman coupling is then turned on adiabatically, as dis- 
cussed above. The initial state being a superposition 
of the two-lowest energy Schrodinger-cat energy eigen- 
states, is now led to oscillate coherently, and collectively, 
between the right- and left-well states. This oscillation is 
analogous to macroscopic quantum tunnelling observed, 
for example, in the Rabi oscillations of superconducting 
qubits 30] . If the Raman coupling is left on for a du- 
ration which is 1/4 of the Bohr period between the even 
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and odd Schrodinger-cat states, the many-body state \^>) 
will evolve as 

\*) ^ S+) J S - ] = 1 ^\N,0) + 1 -^\0,N) (15) 

The Raman beams can be then switched off - with the 
transformation of the state into a momentum-space or 
internal-state-space Schrodinger-cat state, as desired - 
for a duration r, during which a relative phase can ac- 
crue between the two distinguishable portions of the 
Schrodinger-cat state, i.e. 

I*> - ^yV.o) + i±V'^|0, AO (16) 

Here A is the difference in energy between the single- 
particle states to which the right- and left-well states are 
connected when the Raman beams are extinguished. For 
instance, if we are following the implementation of a mea- 
surement of the clock frequency between internal states, 
A = S. 

Finally, this state must be analyzed. In the scheme of 
Bollinger et al. a conventional it/2 pulse is applied, 
and a precise measurement is made of the difference be- 
tween the number of atoms in the \A) state and the \B) 
state - one must establish whether this number differ- 
ence is even or odd, thus requiring single-atom precision 
in the number counting. 

Alternately, one may apply a second "magic beamsplit- 
ter" pulse as before, thus preparing a state 

I - 1 \N, 0) + ii±^ 1\0, N) 

.(17) 

Measurements on this state would detect all atoms in ei- 
ther the right- or left- well states (say, internal states \A) 
or l-B)), with probabilities which vary periodically with 
the free-evolution time t with frequency NA which is N 
times higher than for a single-particle Ramsey-type mea- 
surement. This is the source of the enhanced Heisenberg- 



level sensitivity. This sensitivity persists even in the pres- 
ence of overall number fluctuations, presuming that one's 
goal is to ascertain whether A = 0. However, the correct 
timing required to yield a perfect "magic beamsplitter" 
will depend on N [15] . Thus, attaining the highest pre- 
cision will still require a highly accurate determination 
of the total atom number AT, though perhaps not nec- 
essarily at the single-atom level. A further unexamined 
issue is how the finite temperature of the atomic sample, 
which prevents the initial state from being a pure 
Bose-Einstein condensate, will impact the precision of 
these Schrodinger-cat-based measurements. 

In conclusion, we have presented an analytic model 
which identifies the many-body ground states of a 
weakly-interacting Bose gas which is "periodically 
dressed" by continuous Raman excitation and confined 
in an harmonic spatial potential. The system is ana- 
lyzed in momentum space, wherein the balance between 
tunnelling and weak interactions dictates whether the 
ground states are uncorrelated product states of single- 
particle wavefunctions, or highly correlated states. It is 
found that interactions in the physical system consid- 
ered here have the opposite effect as in the dual sit- 
uation of a position-space double-well potential, that 
is, repulsive interactions will lead to Schrodinger-cat 
states while attractive interactions will lead to number- 
squeezed states with equal numbers of atoms in each well. 
Further, we show that experimental parameters can be 
used to dynamically tune the interaction strength and 
tunnelling rates. This degree of control can be used to 
generate maximally-entangled states directly suitable for 
Heisenberg-limited metrology and interferometry. 
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